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ABSTRACT 



Aims. We search for and characterize substructures in the projected distribution of galaxies observed in the wide field CCD images 
of the 77 nearby clusters of the Wide-field Nearby Galaxy-cluster Survey (WINGS). This sample is complete in X-ray flux in the 
redshift range 0.04 < z < 0.07. 

Methods. We search for substructures in WINGS clusters with DEDICA, an adaptive-kernel procedure. We test the procedure on 
Monte-Carlo simulations of the observed frames and determine the reliability for the detected structures. 

Results. DEDICA identifies at least one reliable structure in the field of 55 clusters. 40 of these clusters have a total of 69 substructures 
at the same redshift of the cluster (redshift estimates of substructures are from color-magnitude diagrams). The fraction of clusters 
with subclusters (73%) is higher than in most studies. The presence of subclusters affects the relative luminosities of the brightest 
cluster galaxies (BCGs). Down to L ~ 10 11 2 L Q > our observed differential distribution of subcluster luminosities is consistent with the 
theoretical prediction of the differential mass function of substructures in cosmological simulations. 

Key words. Galaxies: clusters: general - Galaxies: kinematics and dynamics 



1. Introduction 

According to the current cosmological paradigm, large struc- 
tures in the Universe form hierarchically. Clusters of galaxies 
are the largest structures that have grown through mergers of 
smaller units and have achieved near dynamical equilibrium. In 
the hierarchical scenario, clusters are a rather young population, 
and we should be able to observe their formation process even at 
rather low redshifts. A signature of such process is the presence 
of cluster substructures. A cluster is said to contain substruc- 
tures (or subclusters) when its surface density is characterized 
by multiple, statistically significant peaks on scales larger than 
the typical galaxy size, with "surface density" being referred to 
the cluster gala xies, the intr a-cluster (IC) gas or the dark matter 
(DM hereafter; lBuotell2002l) . 

Studying cluster substructure therefore allows us to investi- 
gate the process by which clusters form, constrain the cosmo- 
logical model of structure for mation, and ultimately test the hi- 
erarchical paradigm itself (e.g. iRichstone et atTl9 92; Mohr et alj 
ll995tlThomas et aljl998l) . In addition, it also allows us to better 
understand the mechanisms affecting galaxy evolution in clus- 
ters, which can be accelerated by the perturbative effects of a 
cluster-subcluster collision and of the tidal field experienced by 
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a group accre ting onto a cluster ( Bekkl IT999I : Ibubinski 1999; 
Gnedin 1999). If clusters are to be used as cosmological tools, 
it is important to calibrate the effects su bstructures have on 
the estimate of their internal properties (e.g . Schindler & Muller 
| 1993[ iPinknev et "all 119961; iRoettiger et alJll998t iBiviano et~ 
2006; lLopes et al.1 l2006l) 7 Finally, detailed analyses of clus- 
ter substructures can be used to constrai n the nature of DM 
dMarkevitch et al.ll2004t Iciowe et alj|2006l) . 

The analysis of cluster substructures can be performed us- 
ing t he projected phase- space distribution of cluster galaxies 
(e.g. iGeller & Beerslll982l) . the surface-brightness distribution 
and t emperature of the X-ray emitting IC gas (e.g. iBriel et all 
1992), or the shear pattern in the background galaxy distribu- 
tion induced by gravitational lensing, t hat directly samples sub - 
structure in the DM component (e.g. Abdelsala m et al.l Fl998). 
None of these tracers of cluster substructure (cluster galaxies, IC 
gas, background galaxies) can be considered optimal. The iden- 
tification of substructures is in fact subject to different biases 
depending on the tracer used. In X-rays projection effects are 
less important than in the optical, but the identification of sub- 
structures is more subject to a z-dependent bias, arising from the 
point spread function of the X -ray telescope and detector (e.g. 
iBohringer & SchueckeiT 2002). Moreover, the different cluster 
components respond in a different way to a cluster-subcluster 
collision. The subcluster IC gas can be ram-pressure braked and 
stripped from the colliding subcluster and lags behind the sub- 
cluster galaxies and DM along the direction of collision (e.g. 
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iRoettiger et all 11997b iBarrena et al.1 l2002t IClowe et alj [2006). 
Hence, it is equally useful to address cluster substructure analy- 
sis in the X-ray and in the optical. 

Traditionally, the first detections of cluster substructures 
were obt ained from the projected s patial distributions of galax- 
ies (e.g. IShane & Wirtanenl 1 1954 lAbell et all 11964 . in com- 
bination, w hen possible, with the distribution of galaxy ve - 
locities (e.g. Ivan den Berghlll960l Il96ll Ide Vaucouleurs|[l96lh . 
Increasingly sophisticated techniques for the detection and 
characterization of cl uster substructures have been develop ed 

19861: iPerea et all Il986allbl 



over the years (see Moles et a l] 



Buoti|2002t iGirardi & Bivianol 12002 " and references therein). 



In many of these techniques substructures are identified as de- 
viations from symmetry in the spatial and/or velocity distri- 
bution of galaxi e s and in the X-ray surf a ce-brightness (e.g . 
West et alJ TT981 iFitchett & Merrittl IT988I: iMohr et all [1993; 



Schuecker et al 



20011) . In other techniques substructures are 



identified as significant peaks in the surface density distribu- 
tion of galaxies or in the X-ray surface brightness, either as 
residuals left after the subtractio n of a smooth, regular model 
representation of the cluster (e.g. Neumann & B ohringer 19971; 
lEttori et aill 19981) . or in a non-parametric way, e.g. by the tech- 
nique of wavelets (e .g. lEscalera et alj [l994; Slezaketal][l994; 
Biviano et al. 1996) and by adaptive-kernel tech niques (e.g. 
Kriessler & Beerslll997tlBardelli et al.lfl998aLl200Tb . 



The performances of several different methods have been 
evalu at ed both usi n g numerical simulati o ns (e. g . iMohr et alJ 
1995 1 ICroneet all Il996t iPinknev etafl Il996t iBuote&Xul 
1997t ICenlll997h IValdarnini et alJll999t iKnebe & Mulleri l2000; 
Biviano et al.l 120061) and also by applying different methods 
to the sa me cluster data-sets and examine the result differ- 
ences (e.g. lEscalera et al.Jll992l Il994t IMohr et alJll995l I" l996; 
Kriessler & Beers 1997; Faddaetal. 1998; Kolokot ronis et alj 
20011: ISchueckeretalJ 120011 lLopes et all 120061) . Generally 



speaking, the sensitivity of substructure detection increases with 
both increasing statistics (e.g. more galaxies or more X-ray pho- 
tons) and increasing dimensionality of the test (e.g. using galaxy 
velocities in addition to their positions, or using X-ray tempera- 
ture in addition to X-ray surface brightness). 

Previous investigations have found very different fractions 
of clusters with substructure in nearby clusters, depending 
on the method and tracer used for substructure detection, on 
the cluster sample, and on the size of sampled clus t er re- 



gions (e.g. iGeller & Beersl 1982; Dressier & Sh ectmai 
etallfTT " 



Moto- 



rs 



1988 



1997 



1995b IG irardi et al.lll997H Kriessler & Been 
Jones & Forman 1999; Solanes et al. 1999; Kolokotronis et al 
200 1 : ISchuecker et alj|200lt iFlin & Krvwultll2006l lLopes et al 
2006). Although the distribution of subcluster masses has not 
been determined observationally, it is known that subclusters 
of ~ 10% the cluster mass are typical, while more massive 
subclusters are less frequent (lEscalera et alJll994t IGirardi et al.1 



119971: I Jones & Formanll 19991) . The situation is probably dif- 
ferent for distant clusters which tend to show massive sub- 
structures more often than nearby clusters clearly suggesting 
hiera rchical growth of clusters was more intense in the past 
(e.g. Gioia et al.l 119991; Ivan Dokkum et all 120001: lHaines et al.l 



2001; Maughan et al. 120031; IHuo et alj|2004t iRosati et al.112004 
Demarco et alj|2005t Ijeltema et al.ll2005l) . 



Additional evidence for the hierarchical formation of clus- 
ters is provided by the analysis of brightest cluster galaxies 
(BCGs hereafter) in substructured clusters. BCGs usually sit 
at the bottom of th e potential well of their host cluster (e.g. 
lAdami et ai1 [l998b). When a BCG is found to be significantly 
displaced from its cluster dynamical center, the cluster displays 



evide nce of substructure (e.g. iBeers et all 1199 It iFerrari et al.l 
1200% From the correlation between cluster and BCG luminosi- 
ties, iLin & Mohrl (120041) conclude that BCGs grow by merg- 
ing as their host clusters grow hierarchically. The related evo- 
lution of BCGs and their host clusters is also suggest ed by the 
alignement of the main cluster and BCG axes (e.g. iBinggelil 
Il982t iDurret et al.lll998l) . Both the BCG and the cluster axes 
are aligned with the surrounding large scale structure dis- 
tribution, where infalling groups come from. These infalling 
groups are finally ident ified as substructures once they enter 
the cluster environment ( Durret et"ailll998t lArnaud et al.l [2000; 
West &Blakesled 120001 [Ferrari et all 120031: iPlionis et alj |2003 



Adam i et al. 120051 Hence, substructure studies really provide 



direct evidence for the hierarchical formation of clusters. 

Concerning the impact of subclustering on global cluster 
properties, it has been found that subclustering leads to over- 
estimating cluster velocity dispersions and virial mas ses (e.g. 
iPerea et alJll99"oT: iBirdll 19951: iMaurogordato et al.ll2000l). but not 
in the general case of small substru ctures (Escalera et al. 1994; 
IGirardi et al.l 119971; IXu et al.l l2000b . During the collision of a 
subcluster with the main cluster, both the X-ray emitting gas 
distribution and its temperature have been found t o be signifi- 
cantly affected (e.g. iMarkevitch & Vikhlini5l200lt IClowe et all 
2006). As a consequence, it has been argued that substruc- 
ture can explain at least part of the scatter in th e scaling rela- 
tions of optical-to-X-ray cluster properties (e.g. Fitchettl ll 988; 



IGirardi et all 1 9961; IBarrena et al.ll2002t lLopes et alJliooi 

As far as the internal properties of cluster galaxies are 
concerned, there is observational evidence that a higher frac- 
tion of cluster galaxies with spectral features characteristic 
of recent or ongoing starburst episodes is located in sub- 
structures or in the regions of cluster-subcluster interactions 
(ICaldwell etafll993lffiraham et alJll996tlBiyiano et al.ll 19971: 
Caldwell & Rosd 119971 iBardelli et all Il998bl iMoss & Whittle! 
2000 : iMiller et al.1 12004 iPoggianti et all l2004t iMillerl 120051; 
Giaci ntucci et al.ll2006l) 



In this paper we search for and characterize substructures 
in the sample of 77 nearby clusters of the Wide-field Nearb y 
Galaxy-cluster Survey (WINGS hereafter. iFasano et al.l l2006). 
This sample is an almost complete sample in X-ray flux in the 
redshift range 0.04 < z < 0.07. We detect substructures from the 
spatial, projected distribution of galaxies in the clus ter fields, us- 
ing th e adaptive-kernel based DEDICA algorithm (Pisani 1993, 
1996). In Sect. [2] we describe our data-set; in Sect. [3] we de- 
scribe the procedure of substructure identification; in Sect. |4] we 
use Monte Carlo simulations in order to tweak our procedure; 
in Sect. [5] we describe the identification of substructures in our 
data-set; in Sect. |6]the catalog of identified substructures is pro- 
vided. In Sect. Q we investigate the properties of the identified 
substructures, and in Sect. [8] we consider the relation between 
the BCGs and the substructures. We provide a summary of our 
work in Sect. [9] 

2. The Data 

WINGS is an all-sky, photometric (multi-band) and spectro- 
scopic survey, whose global goal is the systematic study of the 
local cosmic variance of the cluster population and of the prop- 
erties of cluster galaxies as a function of cluster properties and 
local environment. 

The WINGS sample consists of 77 clusters selected from 
three X-ray flux limited samples compiled from ROSAT All- 
Sky Survey data, with constraints just on the redshift (0.04 < z < 
0.07) and distance from the galactic plane {\b\ >20 deg). The core 
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of the project consists of wide-field optical imaging of the se- 
lected clusters in the B and V bands. The imaging data were col- 
lected using the WFC@INT (La Palma) and the WFI@MPG (La 
Silla) in the northern and southern hemispheres, respectively. 

The observation strategy of the survey favors the uniformity 
of photometric depth inside the different CCDs, rather than com- 
plete coverage of the fields that would require dithering. Thus, 
the gaps in the WINGS optical imaging correspond to the phys- 
ical gaps between the different CCDs of the mosaics. 

During the data reduction process, we give particular care to 
sky subtraction (also in presence of crowded fields including big 
halo galaxies and/or very bright stars), image cleaning (spikes 
and bad pixels) and star/galaxy classification (obtained with both 
automatic and int eractive tools) . 

According to iFasano et ail d2006l) and IVarela et all d2007l) . 
the overall quality of the data reported in the WINGS photomet- 
ric catalogs can be summarized as follows: (i) the astrometric 
errors for extended objects have r.m.s. ~0.2 arcsec; (ii) the av- 
erage limiting magnitude is ~24.0, ranging from 23.0 to 25.0; 
(iii) the completeness of the catalogs is achieved (on average) up 
to V ~22.Q; (iv) the total (systematic plus random) photometric 
r.m.s. errors, derived from both internal and external compar- 
isons, vary from ~0.02 mag, for bright objects, up to ~0.2 mag, 
for objects close to the detection limit. 

3. The DEDICA Procedure 

We base our search for substructures in WINGS clusters on the 
DEDICA procedure dPisanill 1 993l 1 1 9961) . This procedure has the 
following advantages: 

1 . DEDICA gives a total description of the clustering pattern, 
in particular the membership probability and significance of 
structures besides geometrical properties; 

2. DEDICA is scale invariant; 

3. DEDICA does not assume any property of the clusters, i.e. 
it is completely non-parametric. In particular it does not re- 
quire particularly rich samples to run effectively. 

The basic nature a nd properties of DEDICA are described in 
Pisanl ( 1 19931 [T996, and references therein). Here we summarize 
the main structure of the algorithm and how we apply it to our 
data sample. The core structure of DEDICA is based on the as- 
sumption that a structure (or a "cluster" in the algorithm jargon) 
corresponds to a local maximum in the density of galaxies. 

We proceed as follows. First we need to estimate the proba- 
bility density function ^(r,) (with i = 1, . . .AO associated with 
the set of N galaxies with coordinates r,. Second, we need to 
find the local maxima in our estimate of T(r,-) in order to iden- 
tify clusters and also to evaluate their significance relatively to 
the noise. Third and finally, we need to estimate the probability 
that a galaxy is a member of the identified clusters. 



3. 1 . The probability density 

DEDICA is a non-parametric method in the sense that it does 
not require any assumption on the probability density function 
that it is aimed to estimate. The only assumptions are that ^(r,) 
must be continuous and at least twice differentiable. 

The function /(r,) is an estimate of *T(r,) and it is built by 
using an adaptive kernel method given by: 



1 N 

Mr) = — ^ K(r u err, r) 

i=l 



where we use the two dimensional Gaussian kernel K(r h <t(, r) 
centered in r, with size cr,-. 

The most valuable feature of DEDICA is the procedure to se- 
lect the values of kernel widths cr, . It is possible to show that the 
optimal choice for cr,, i.e. with asymptotically minimum vari- 
ance and null bias, is obtained by minimizing the distance be- 
tween our estimate /(r,) and ( P(r,). This distance can be eval- 
uated by a particular function called the integrated square error 
ISE(f) given by: 



ISE(f) 



[<F(r)-/(r)] 2 rfr 



(2) 



Once the minimum ISE( f) is reached we have obtained the 
DEDICA estimate of the density as in EqQ] 



3.2. Cluster Identification 

The second step of DEDICA consists in the identification of the 
local maxima in /t a (r). The positions of the peaks in the density 
function fka(r) are found as the solutions of the iterative equa- 
tion: 



r m + i 



r m + a 



(3) 



where a is a scale factor set according to optimal convergence 
requirements. The limit R of the sequence r,„ defined in Eqf3] 
depends on the starting position r m= i. 



lim r„ 

m— >+oo 



R(r m= i) 



(4) 



(1) 



We run the sequence in Eqf5]at each data position r,. We label 
each data point with the limit R, = R(r m= i = r,). These limits 
R, are the position of the peak to which the i - th galaxy belong. 
In the case that all the galaxies are members of a unique cluster, 
all the labels R, are the same. At the other extreme each galaxy 
is a one-member cluster and all R, have different values. All the 
members of a given cluster belong to the same peak in //t fl (r) 
and have the same R,. We identify cluster members by listing 
galaxies having the same values of R. We end up with v different 
clusters each with (p — 1, . . ., v) members. 

In order to maintain a coherent notation, we identify with 
the label fi — the «o isolated galaxies considered a system of 
background galaxies. We have: hq - N - YL=\ n n- 

3.3. Cluster Significance and of Membership Probability 

The statistical significance (fi — 1, ... , v) of each cluster is 
based on the assumption that the presence of the /i - th cluster 
causes an increase in the local probability density as well as in 
the sample likelihood = n,-[/* fl (r,-)] relatively to the value 
that one would have if the members of the fi - th cluster were all 
isolated, i.e. belonging to the background. 

A large value in the ratio Ln/L u characteri zes the most im- 
portant clusters. According to Materne] d 1979b it is possible to 
estimate the significance of each cluster by using the likelihood 
ratio test. In other words 2\n{L^ I L^) is distributed as ax 2 vari- 
able with v - 1 degrees of freedom. Therefore, once we compute 
the value of x 1 for each cluster (x\), we can also compute the 
significance S ^ of the cluster. 

Here we assume that the contribution to the global density 
field /i fl (r,) of the fi - th cluster is F^r,-). The ratio between the 
value of F^r,) and the total local density /A fl (r,) can be used to 
estimate the membership probability of each galaxy relatively to 
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the identified clusters. This criterion also allows us to estimate 
the probability that a galaxy is isolated. 

At the end of the DEDICA procedure we are left with a) 
a catalog of galaxies each with information on position, mem- 
bership, local density and size of the Gaussian kernel, b) a cat- 
alog of structures with information on position, richness, the 
^parameter, and peak density. For each cluster we also com- 
pute from the coordinate variance matrix the cluster major axis, 
ellipticity and position angle. 

4. Tweaking the Algorithm with Simulations 

In this section we describe our analysis of the performance of 
DEDICA and the guidelines we obtain for the interpretation of 
the clustering analysis of our observations. 

We build simulated fields containing a cluster with and with- 
out subclusters. The simulated fields have the same geometry 
of the WFC field and are populated with the typical number of 
objects we will analyze. For simplicity we consider only WFC 
fields. Because DEDICA is scale-free, a different sampling of 
the same field of view has no consequence on our analysis. 

In the next section we limit our analysis to Myji,,, < -16. 
At the median redshift of the WINGS cluster, z 0.05, this 
absolute magnitude limit corresponds to an apparent magnitude 
Vum -21. Within this magnitude limit the representative number 
of galaxies in our frames is Ntot— 900. 

We then consider N tot - N mem + Nbkg, with N mem \ht number 
of cluster members and iV^the number of field - or background 
- galaxies. We set Nbkg- 670, close to the average number of 
background galaxies we ex pect in our frame ba sed on typical 
observ ed fields counts, e.g. lBerta et alj d2006l) or lArnouts et alj 
d 1997b . With this choice, we have N mem = 230. 

We distribute uniformly at random A^MgObjects. We dis- 
tribute at random the remaining N mem - 230 objects in one or 
more overdensities depending on the test we p erform. We p opu- 
late overdensities according to a King profile dKing| [l962) with 
a core radius R core = 90 kpc, representative of our clusters. We 
then scale R core with the number of members of the substructure, 
Ns ■ We use 



where Nc is the number of objects in the cluster with N mem — 
Ns + Nc- T h is scal ing of R core with cluster richness is from 
AdamTetaD £[998a) assuming direct proportionality between 
cluster richness and luminosity (e.g. lPopesso et al]|2006l) . 

As far as the relative richnesses of the cluster and subcluster 
are concerned, we consider the following richness ratios r cs = 
Nc/Ns = 1,2,4,8. With these richness ratios, the number of 
objects in the cluster are Nc = 115, 153, 184,204, and those in 
subclusters are Ns = 115, 77, 46, 26 respectively. 

In a first set of simulated fields we place the substructure at 
2731 pixels (15 arcmin) from the main cluster so that they do not 
overlap. In a second set of simulations, we place main cluster and 
substructure at shorter distances, 683 and 1366 pixels, in order to 
investigate the ability of DEDICA to resolve structures. At each 
of these shorter distances we build simulations with both r cs = 1 
and 2. 

For each richness ratio and/or distance between cluster and 
subcluster we produce 16 simulations with different realizations 
of the random positions of the data points representing galaxies. 

In order to minimize the effect of the borders on the detection 
of structures we add to the simulation a "frame" of 1000 pixel. 




simulation 



Fig. 1. Fraction of recovered members of each substructure for 
different r cs . The solid line connects substructures with r cs = 2 
and 4 

We fill this frame with a grid of data points at the same density 
as the average density of the field. 

The first result we obtain from the runs of DEDICA on the 
simulations with varying richness ratio is the positive rate at 
which we detect real structures. We find that we always recover 
both cluster and substructure even when the substructure only 
contains Ns = 1/8 Nc objects, i.e. 26 objects (on top of the uni- 
form background). In other words, if there is a real structure 
DEDICA finds it. 

We also check how many original members the procedure 
assigns to structures it recovers. The results are summarized in 
Fig-E In the diagram, the fraction of recovered members of each 
substructure is represented by the values of its r cs . The solid line 
connects substructures with r cs = 2 and 4. 

From Fig. [TJ it is clear that our procedure recovers a large 
fraction of members, almost irrespective of the richness of the 
original structure. It is also interesting to note that the fluctu- 
ations identified as substructures are located very close to the 
center of the corresponding simulated substructures. In almost 
all cases the distance between original and detected substructure 
is significantly shorter than the mean inter-particle distance. 

The second important result we obtain from the simulations 
is the false positive rate, i.e. the fraction of noise fluctuations 
that are as significant as the fluctuations corresponding to real 
structures. 

First of all we need to define an operative measure of the 
reliability of the detected structures. In fact DEDICA provides a 
default value 5^ (/y = 1, . . . , v) of the significance (see Sect. l3.3b - 
However, S M has a relatively small dynamical range, in particular 
for highly significant clusters. 

Density or richness both allow a reasonable "ranking" of 
structures. However, both large low-density noise fluctuations 
(often built up from more than one noise fluctuation) and very 
high density fluctuations produced by few very close data points 
could be mistakenly ranked as highly significant structures ac- 
cording to, respectively, richness and density criteria. 
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20 40 60 



xi 

Fig. 2. x 2 s of simulated noise fluctuations (solid line). Labels are 
the r cs of simulated structures at the abscissa corresponding to 
their x\ an d at arbitrary ordinates. 



We therefore prefer to use the parameter^ which stands at 
the base of the estimate of 5^ and which is naturally provided by 
DEDICA. The main characteristic of^is that it depends both on 
the density of a cluster relative to the background and on its rich- 
ness. Using x\ we classify correctly significantly more structures 
than with either density or richness alone. 

In Fig. |2] we plot the distribution of^of noise fluctuations 
(solid line). In the same plot we also mark the r f 5 of real struc- 
tures as detected by our procedure. We use labels indicating 
r„and place them at the abscissa corresponding to their x\ an d 
at arbitrary ordinates. 

Fig.|2]shows that the structures detected with r cs - 1, 2 are al- 
ways distinguishable from noise fluctuations. Substructures with 
r cs - 4 or higher, although correctly detected, have x\ values that 
are close to or lower than the level of noise. 

With the second set of simulations, we test the minimum dis- 
tance at which cluster and subcluster can still be identified as 
separate entities. We place cluster and substructure (r cs - 1, 2) at 
distances d cs = 683 and 1366 pixel. These distances are 1/4 and 
1/2 respectively of the distance between cluster and substructure 
in the first set of simulations. Again we produce 16 simulations 
for each of the 4 cases. 

We find that at d cs = 1366 pixel cluster and substructure are 
always correctly identified. At the shorter distance d cs = 683 
pixel, DEDICA merges cluster and substructure in 1 out of 16 
cases for r cs - 1 and in 8 out of 16 cases for r cs - 2. With our 
density profile, d cs = 683 pixel corresponds to d cs ^ R c + R s 
with R c , R s the radii of the main cluster and of the subcluster 
respectively. 

In order to verify the results we obtain for 900 data points 
we produce more simulations with N tot - 450, 600 and 1200. In 
all these simulations Rc and Rs are the same as in the set with 
N to t= 900. We vary A^and N me ,„so that N mem / Nbkgis the same 
as in the case N tot - 900. 



X 




0.5 1 1.5 2 2.5 
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Fig. 3. Small symbols correspond tox\ as a function of the num- 
ber of members of noise fluctuations. Crosses, circles, dots and 
triangles are x\ for the noise fluctuations of the simulations with 
Ntot- 450, 600, 900, and 1200 respectively. Large symbols are 
XgOf simulated clusters and subclusters with r cs - 1. Horizontal 
lines mark the levels of x\ threshold- 



These simulations confirm the results we obtain in the 
case Ntot- 900, and allow us to set a detection threshold, 
X~s threshold^ tot ^ ^ or significant fluctuations in the analysis of real 
clusters. 

We summarize the behavior of the noise fluctuations in our 
simulations in Fig. [3] In this figure, the small symbols corre- 
spond to x\ as a function of the number of members of noise 
fluctuations. In particular, crosses, circles, dots and triangles are 
^for the noise fluctuations of the simulations with Ntot— 450, 
600, 900, and 1200 respectively. 

The larger symbols are the x\ °f the fluctuations correspond- 
ing to simulated clusters and subclusters of equal richness (r cs - 
1). 

The 4 horizontal lines mark the level of x\ threshold' tne 
average x\ °f the 3 most significant noise fluctuations in each of 
the 4 groups of simulations with N t0 ,= 450, 600, 900, and 1200. 

The expected increase ofx s threshold™^ 1 NtoA$ evident. 

We note that the only significant difference with these find- 
ings we obtain from the simulations with r cs - 2 is that Jtjof 
simulated clusters and subclusters is closer to^ thre ^(but still 
higher). 

We m xi,hreshold™ th ^'»' and ° btain 

^(xhhreshoJ = 1 /2-55 log(AW + 0.394 (5) 

in good agreement with the expected behavior of the poissonian 
fluctuations. 

As a final test we verify that infra-chip gaps do not have a 
dramatic impact on the detection of structures in the cases r cs - 1 
and 2. We place a 50 pixel wide gap where it has the maximum 
impact, i.e. where the kernel size is shortest. Even if the infra- 
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chip gap cuts through the center of the structures, DEDICA is 
able to identify these structures correctly. 

We summarize here the main results of our tests on simulated 
clusters with substructures: 

- DEDICA successfully detects even the poorest structures 
above a uniform poissonian noise background. 

- DEDICA recovers a large fraction (typically > 3/4) of the 
real members of a substructure, almost irrespective of the 
richness of the structure. 

- DEDICA is able to distinguish between noise fluctuations 
and true structures only if these structures are rich enough. 
In the case of our simulations, structures have to be richer 
than l/4th of the main structure. 

- DEDICA is able to separate neighboring structures provided 
they do not overlap. 

- infra-chip gaps do not threaten the detection of structures 
that are rich enough to be reliably detected. 

- the x\ threshold we use to identify significant structures is 
a function of the total number of points and can be scaled 
within the whole range of numbers of galaxies observed 
within our fields. 




In the next section we apply these results to the real WINGS 
clusters. 



5. Substructure detection in WINGS clusters 

We apply our clustering procedure to the 77 clusters of the 
WINGS sample. The photometric catalog of each cluster is deep, 
reaching a completeness magnitude V comp i e , e < 22. The num- 
ber of galaxies is correspondingly large, from N ga i - 3, 000 to 
Ngal =* 10,000. 

The large number of bright background galaxies (faint ap- 
parent magnitudes) dilutes the clustering signal of local WINGS 
clusters. We perform test runs of the procedure on several clus- 
ters with magnitude cuts brighter than V complete- Based on these 
tests, we decide to cut galaxy catalogs to the absolute magnitude 
threshold My = -16.0. With this choice a) we maximize the 
signal-to-noise ratio of the detected subclusters and b) we still 
have enough galaxies for a stable identification of the system. At 
the median redshift of WINGS clusters, z — 0.0535, our absolute 
magnitude cut corresponds to an apparent magnitude V - 21.2. 

This apparent magnitude also approximately corresponds to 
the magnitude where the contrast of our typical cluster relative 
to the field is maximum (this estimate is based on the average 



cluste r luminosity function of dYagi et al.l2002| De Propris et alj 
120031) and on the galaxy counts of dBerta et all2006h ). 

The number of galaxies that are brighter than the threshold 
M v = -16.0 is in the range 600 < N, ol < 1200 for a large 
fraction of clusters observed with either WFC@INT or with 
WFI@ES02.2. 

In order to proceed with the identification of significant 
structures within WINGS clusters, we need to verify that our 
simulations are sufficiently representative of the real cases. 
In practice we need to compare the observed distributions of 
Rvalues of noise fluctuations with the corresponding simulated 
distributions. In the observations it is impossible to identify indi- 
vidual fluctuations as noise. In order to have an idea of the distri- 
butions of ^of noise fluctuations we consider that our fields are 
centered on real clusters. As a consequence, on average, fluctua- 
tions in the center of the frames are more likely to correspond to 
real systems than those at the borders. 



Fig. 4. ^distributions for border (thick solid histogram) and 
central (thick dashed histogram) observed fluctuations. The thin 
solid line is the normalized distribution of^of the noise fluctu- 
ations in our simulations 



We therefore consider separately the fluctuations within the 
central regions of the frames and all other fluctuations (borders). 
We define the central regions as the central 10% of WFC and 
WFI areas. We plot in Fig. [4] the two distributions. The thick 
solid histogram is for the border and the thick dashed histogram 
for the center of the frames. The difference between "noise" and 
"signal" is clear. In the same figure we also plot the normal- 
ized distribution of^of the noise fluctuations in our simulations 
(thin solid line). The distributions of^of the observed and sim- 
ulated fluctuations are in reasonable agreement considering a) 
the simple model used for the simulations and that b) in the ob- 
servations we can not exclude real low-;f ^ structures among noise 
fluctuations. We conclude that for our clusters we can adopt the 
same reliability threshold \ 2 S ,hresiwid we determine from our sim- 



ulations (Eq.|5]l. 



6. The Catalog of Substructures 

We detect at least one significant structure in 55 (71%) clus- 
ters. We find that 12 clusters (16%) have no structure above the 
threshold (undetected). In the case of another 10 (13%) clusters 
we find significant structures only at the border of the field of 
view. In absence of a detection in the center of the frame, we con- 
sider these border structures unrelated to the target cluster. We 
also verify that in the Color-Magnitude Diagram (CMD) these 
border structures are redder than expected given the redshift of 
the target cluster. We consider also these 10 clusters undetected. 

Here we list the 22 undetected clusters: A0133, A0548b, 
A0780, A1644, A1668, A1983, A2271, A2382, A2589, A2626, 
A2717, A3164, A3395, A3490, A3497, A3528a, A3556, A3560, 
A3809, A4059, RX1022, Z1261. 

We note that undetected clusters are real physical systems ac- 
cording to their x-ray selection. From an operative point of view, 
the fact that these clusters are not detected by DEDICA is the 
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Fig. 5. Isodensity contours (logarithmically spaced) of the Abell 
85 field. The title lists the coordinates of the center. The orien- 
tation is East to the left, North to the top. Galaxies belonging to 
the systems detected by DEDICA are shown as dots of different 
colors. Black, light green, blue, red, magenta, dark green are for 
the main system and the subsequent substructures ordered as in 
Table lA.ll Large symbols are for galaxies with My < -17.0 that 
lie where local densities are higher than the median local density 
of the structure the galaxy belongs to. Open symbols mark the 
positions of the first- and second-ranked cluster galaxies, BCG1 
and BCG2 respectively. Similar plots for the 55 analysed clusters 
are available in the electronic version of this Journal. 



result of the division into too many structures of the total avail- 
able clustering signal in the field (or of a too large fraction of 
the clustering signal going into border structures). Several phys- 
ical situations could be at the origin of missed detections. One 
possibility is an excess of physical substructures of comparable 
richness. Another possibility is that these clusters are embedded 
in regions of the large scale structure that are highly clustered. 

We do not try to recover these structures because they can 
not be prominent enough. Since our analysis is bidimensional, 
we can only detect and use confidently the most prominent struc- 
tures. Redshifts are needed for a more detailed analysis of cluster 
substructures. 

We list the 55 clusters with significant structures in Table 
I A. 1 1 We give, for each substructure: (1) the name of the parent 
cluster; (2) the classification of the structure as main (M), sub- 
cluster (S), or background (B) together with their order number; 
(3) right ascension (J2000), and (4) declination (J2000) in deci- 
mal degrees of the DEDICA peak; the parameters of the ellipse 
we obtain from the variance matrix of the coordinates of galaxies 
in the substructure, i.e. (5) major axis in arcminutes, (6) elliptic- 
ity, and (7) position angle in degrees; (8) luminosity (see the next 
section); (9)x 2 s ■ 

We make available contour plots of the number density fields 
of all clusters in Fig.|6]of the electronic version of this Journal. 
In Fig.|5]we show an example of these plots. Isodensity contours 
are drawn at ten logarithmic intervals. Galaxies belonging to the 



systems detected by DEDICA are shown as dots of different col- 
ors. We use large symbols for brighter galaxies (My < -17.0) 
that lie where local densities are higher than the median local 
density of the structure the galaxy belongs to. We also mark 
with open symbols the positions of the first- and second-ranked 
cluster galaxies, BCG1 and BCG2 respectively. Color coding is 
black, light green, blue, red, magenta, dark green for the main 
system and the subsequent substructures ordered as in Table lA. fl 
We describe and analyze in detail our catalog in the next sec- 
tion. 

7. Properties of substructures 

The first problem we face in order to study the statistical and 
physical properties of substructures is to determine their asso- 
ciation with the main structure. In fact, the main structure itself 
has to be identified among the structures detected by DEDICA 
in each frame. 

In most cases it is easy to identify the main structure of a 
cluster since it is located at the center of the frame and it has 
a high^. In two cases (AO 168 and A 1736) the choice of the 
main structure is complicated because there are several similar 
structures near the center of the frame. In these cases we select 
the main structure for its highest x\ ■ 

At this point we limit our analysis to members of the struc- 
ture that a) have an absolute magnitude My < -17 (corrected for 
Galactic absorption) and that b) are in the upper half of the distri- 
bution of DEDICA-defined local galaxy densities of the system 
they belong to. The galaxy density threshold we apply allows us 
to separate adjacent structures whose definition becomes more 
uncertain at lower galaxy density levels. The magnitude cut in- 
creases the relative weight of the galaxies we use to evaluate the 
nature of structures in the CMD. 

After having identified the main structure, we need to deter- 
mine which structures in the field of view of a given cluster have 
to be considered background structures. We consider a structure 
a physical substructure (or subcluster) if its color-magnitude re- 
lation (CMR hereafter) is identical, within the errors, to the CMR 
of the main structure. 

As a first step we define the color-magnitude relation (CMR) 
of the "whole cluster", i.e. of galaxies in the main structure to- 
gether with all other galaxies not assigned to any structure by 
DEDICA. We compu te the (B - V) CM R of the Coma cluster 
from published data (Adam i"et alj 12006b . Then we keep fixed 
the slope of the linear CMR of Coma and shift it to the mean 
redshift of the cluster. 

In order to determine that the main structure and a substruc- 
ture are at the same redshift, we evaluate the fraction of back- 
ground (red) galaxies, fi, g , that each structure h as in the CMD . 
If these fractions are identical within the errors dGehrels| [l986). 
we consider the two structures to be at the same redshift. 

In practice we determine /& g by assigning to the background 
those galaxies of a structure that are redder than a line parallel to 
the CMR and vertically shifted (i.e. redwards) by 2.33 times the 
root-mean square of the colors of galaxies in the CMR. We note 
that the probability that a random variable is greater than 2.33 in 
a Gaussian distribution is only 1%. 

The result of the selection of main structures and substruc- 
tures is the following: 40 clusters have a total of 69 substruc- 
tures at the same redshift as the main structure, only 15 clus- 
ters are left without substructures. A total of 35 systems are 
found in the backgro und. Considering a) the number density o f 
poor-to-rich clusters dMazure et alJl996tlZabludoff et al.ll 19931) . 
b) the average luminosity function of clusters dYagi et alj |2002; 
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Fig. 7. Cumulative distributions of the two different indicators of 
subclustering: left panel N sl ,b, right panel fuub- 



Fig. 8. Observed differential distribution of subcluster lumi- 
nosities (histogram) and theoretical model (arbitrary scaling; 
iDe Lucia 'et alj2004 . 



iDe Propris et al.ll2.Q03b . c) the total area covered by the 55 clus- 
ter fields, and d) the limiting apparent magnitude corresponding 
to our absolute magnitude threshold My = -16.0, we expect to 
find ~ 0.5 + 0.2 background systems per cluster field, 28 + 1 1 in 
total. This estimate is consistent with the 35 background systems 
we find. 

The fraction of clusters with subclusters (73%) is higher than 
gen erally found in previous investigations (typically ~ 30%, see, 
e.g., Girardi & Bivianol2002tlFlin & Krvwultl2006tlLopes et al.l 
2006, and references therein). Even if we count all undetected 
clusters as clusters without substructures, this fraction only de- 
creases to 52% (40/77). It is however acknowledged that the 
fraction of substructured clusters depends, among other factors, 
on the algorithm used to detect substructu res, on the quality 
and de pth of the galaxy catalog. For example Kolokotronis et al.l 
(2001) using optical and X-ray dat a find that the fraction of clus- 
ters with substructures is > 45%, Burgettet al. (2004) using a 
battery of tests detect substructures in 84% of the 25 clusters of 
their sample. 

Having established the "global" fraction of substructured 
clusters, we now investigate the degree of subclustering of in- 
dividual clusters, i.e. the distribution of the number of substruc- 
tures N s „b we find in our sample. 

We find 15 (27%) clusters without substructures; 22 (40%) 
clusters wiHhN sub = 1; 10(18%) clusters withA^ = 2;6(11%) 
clusters with N su t = 3; and 2 (3%) clusters with N su b = 4. We 
plot in the left panel of Fig.[7]the integral distribution of N su i,. 

The distribution of the level of subclustering does not change 
when we measure it as the fractional luminosity of subclusters, 
fisub, relative to the luminosity of the whole cluster (see Fig. [7] 
right panel). The lumino sities we est i mate are background cor- 
rected using the counts of lBerta et al J (12006). We use the ellipses 
output from DEDICA (see previous section) as a measure of the 
area of subclusters. 

We find that N su b and fi su j, are clearly correlated according 
to the Spearman rank-correlation test. 

We now consider the distribution of subcluster luminosities 
and plot the corresponding histogram in Fig. [8] In the same fig- 
ure we also plot with arbitrary scaling the power-law oc L l . This 
relation is the prediction for the differential m ass function of 
substr uctures in the cosmological simulations o f De Lucia et al.l 
(2004). 



Our observations are consiste nt to within the uncert ainties 
with the theoretical prediction of IDe Lucia et a D (12004 down 
to L ~ 10" 2 L Q . The disagreement at lower luminosity is ex- 
pected since: a) below this limit galaxy-sized halos become im- 
portant among the simulated substructures, and b) only above 
this limit we expect our catalog to be complete. In fact only 
subclusters with luminosities brighter than L = 10 1L2 L Q have 
always richnesses that are > 1/3 of the main structure. This 
richness limit approximately corresponds to the completeness 
limit of DEDICA detections according to our simulations (see 
Sect.!). 



8. Brightest Cluster Galaxies 

Here we investigate the relation between BCGs and cluster struc- 
tures. 

We find that, on average, BCGls are located close to the den- 
sity peak of the mai n structures. In pro jection on the sky, the bi- 
weight average (see Beers et al. 1990)" distance of BCGls from 
the peak of the main system is 72 ±11 kpc. If we only consider 
the 44 BCGls that are on the CMR and are assigned to main 
systems by DEDICA, the average distance decreases to 56 + 8 
kpc. The fact that BCGls are close to the center of the system 
is consiste nt with current theoretical view o n the formation of 
BCGs (e.g. lDubinskilll998l:lNipoti et al.ll2004l) . 

BCG2s are more distant than BCGls from the peak of the 
main system: the biweight average distance is 345 + 47 kpc. If 
we only consider the 26 BCG2s that are on the CMR and are 
assigned to main systems by DEDICA, the average distance de- 
creases to 161 +34 kpc. 

Projected distances of BCG2s from density peaks remain 
larger than those of BCGls even when we consider the density 
peak of the structure or substructure they belong to. In Fig. [9] 
we plot the cumulative distributions of the distances of BCGls 
(solid line) and BCG2s (dashed line) from the density peak of 
their systems. The distributions are different at the > 99.99% 
level according to a Kolmogorov-Smirnov test (KS-test). 

Now we turn to luminosities and find that the magnitude dif- 
ference between BCGls and BCG2s, AM\i, is larger in clus- 
ters without substructures than in clusters with substructures. 
In Fig. [10] we plot the cumulative distributions of AMn for 
clusters with (dashed line) and without (solid line) subclusters. 
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Fig. 9. Cumulative distributions of distances of BCG1 (solid 
line) and BCG2 (dashed line) from the density peak of their sys- 
tem. 
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Fig. 10. Cumulative distributions of the magnitude difference be- 
tween BCG1 and BCG2 in clusters with (dashed line) and with- 
out subclusters (solid line). 

The two distributions are different acc ording to a KS-test a t the 
99.1% confidence level. We note that iLin & Mohrl d2004 find 
that AM] 2 is independent of cluster properties. These authors 
however do not consider subclustering. 

In order to determine whether the higher values of AM12 in 
clusters without subclusters are due to an increased luminosity 
of the BCG1 (Li) or to a decreased luminosity of the BCG2 
(L2), we consider the luminosity of the 10'' 1 brightest galaxy 
(Z40) as a reference. The biweight average luminosity ratios are 

< L1/L10 >= 8.6 +1.0 and < L 2 /L 10 >= 3.3 ± 0.3 in clus- 
ters without substructures, and < L\/L\q >= 7.1 ± 0.4 and 

< L2/L10 >= 3.4 + 0.2 in clusters with substructures. We then 
conclude that the AMi 2 -effect is caused by a brightening of the 
BCG1 relative to the BCG2 in clusters without substructures. 

The fact that AM12 is higher in clusters without substruc- 
tures can be interpreted, at least qualitatively, in the framework 
of the hierarchical scenario of structure evolution. Clusters with- 
out substructures are likely to be evolved after several merger 
phases. Their BCGls have already had time to accrete many 
galaxies, in particular the more massive ones, which slow down 
and sink to the cluster center as the result of dynamical friction. 



Some of these galaxies may even h ave been BCGs of the merg- 
ing structures. The simulations by De Lucia & Blaizot (2006) 
show that the BCGls continue to increase their mass via can- 
nibalism even at recent times, and that there is a large vari- 
ance in the mass accretion history of BCGls from cluster to 
cluster. The result of such a cannibalism process is an increase 
of the BCG1 luminosity with respect to other cluster galaxies, 
and in extreme cases may lead to the formation of fossil groups 
dKhosroshahi et al.ll2006h . 

However, according to these simulations, only 15% of all 
BCGls have accreted > 30% of their mass over the last 2 Gyr, 
while another 15% have accreted <3% of their mass over the 
same period. Our results indicate that about 60% of the BCGls 
are more than 1 magnitude brighter than the corresponding 
BCG2s. Given the size and generality of the luminosity dif- 
ferences it would seem that cannibalism alone, even if present 
along the merging history of a given cluster, cannot account for 
it. Most of the BCGls should have then been assembled in early 
times, a s pointed out in the downsizing scenario for galaxy for- 
mation dCowie et alJI 19961) and entered that merging history al- 
ready with luminosity not far form the present one. 

9. Summary 

In this paper we search for and characterize cluster substructures, 
or subclusters, in th e sample of 77 nearby clusters of the WINGS 
dFasano et alj |2006). This sample is an almost complete sample 
in X-ray flux in the redshift range 0.04 < z < 0.07. 

We detect substructures in the spatial projecte d distribution 
of ga laxies in the cluster fields using DEDICA dPisanilll993l 
1996) an adaptive-kernel technique. DEDICA has the following 
advantages for our study of WINGS clusters: 

a) DEDICA gives a total description of the clustering pat- 
tern, in particular membership probability and significance of 
structures besides geometrical properties. 

b) DEDICA is scale invariant 

c) DEDICA does not assume any property of the clusters, i.e. 
it is completely non-parametric. In particular it does not require 
particularly rich samples to run effectively. 

In order to test DEDICA and to set guidelines for the in- 
terpretation of the results of the application of DEDICA to our 
observations we run DEDICA on several sets of simulated fields 
containing a cluster with and without subclusters. 

We find that: a) DEDICA always identifies both cluster and 
subcluster even when the substructure richness ratio cluster- 
to-subcluster is r cs = 8, b) DEDICA recovers a large fraction 
of members, almost irrespective of the richness of the original 
structure (k, 70% in most cases), c) structures with richness ra- 
tios r cs £ 3 are always distinguishable from noise fluctuations of 
the poissonian simulated field. 

These simulations also allow us to define a threshold that we 
use to identify significant structures in the observed fields. 

We apply our clustering procedure to the 77 clusters of the 
WINGS sample. We cut galaxy catalogs to the absolute magni- 
tude threshold My = -16.0 in order to maximize the signal-to- 
noise ratio of the detected subclusters. 

We detect at least one significant structure in 55 (71%) clus- 
ter fields. We find that 12 clusters (16%) have no structure above 
the threshold (undetected). In the remaining 10 (13%) clusters 
we find significant structures only at the border of the field of 
view. In absence of a detection in the center of the frame, we 
consider these border structures unrelated to the target cluster. 
We also verify that in the CMD these border structures are redder 
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than expected given the redshift of the target cluster. We consider 
also these clusters undetected. 

We provide the coordinates of all substructures in the 55 
clusters together with their main properties. 

Using the CMR of the early-type cluster galaxies we sep- 
arate "true" subclusters from unrelated background structures. 
We find that 40 clusters out of 55 (73%) have a total of 69 sub- 
structures with 15 clusters left without substructures. 

The fraction of clusters with subclusters (73%) we identify is 
higher tha n most previously publi shed values (typically ~ 30%, 
see, e.g.. iGirardi & Bivianol l2002. and references therein). It is 
however acknowledged that the fraction of substructured clus- 
ters depends, among other factors, on the algorithm used to de- 
tect substructures, on the quality and depth o f the galaxy catalog 
dKolokotronis et alj|200ll;feurgett et aLlliooH) . 

Another important result of our analysis is the distribution of 
subcluster luminosities. In the luminosity range where our sub- 
structure detection is complete (L > 1O 1L2 L ), we find that the 
distribution of subcluster luminosities is in agreement with the 
power-law oc L _1 predicted for the differential mass function of 
substr uctures in the cosmological simulations o f lDe Lucia et alj 
(120041) . 

Finally, we investigate the relation between BCGs and clus- 
ter structures. 

We find that, on average, BCGls are located close to the den- 
sity peak of the main structures. In projection on the sky, the bi- 
weight average distance of BCGls from the peak of the main 
system is 72 + 1 1 kpc. BCG2s are significantly more distant than 
BCGls from the peak of the main system (345 + 47 kpc). 

The fact that BCGls are close to the center of the system 
is cons istent with curr ent theoretical view on the formation of 
BCGs dDubinskilll998l) . 

A more surprising result is that the magnitude difference be- 
tween BCGls and BCG2s, AM12, is significantly larger in clus- 
ters without substructures than in clusters with substructures. 
This fact may be interpreted in the f ramework of the hierarchic al 
scenario of structure evolution (e.g. lDe Lucia & Blaizot 2006). 
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Appendix A: The catalog of substructures 

We provide here the catalog of substructures. In Table I A. ll we give, for each substructure: (1) the name of the parent cluster; (2) the 
classification of the structure as main (M), subcluster (S), or background (B) together with their order number; (3) right ascension 
(J2000), and (4) declination (J2000) in decimal degrees of the DEDICA peak; the parameters of the ellipse we obtain from the 
variance matrix of the coordinates of galaxies in the substructure, i.e. (5) major axis in arcminutes, (6) ellipticity, and (7) position 
angle in degrees; (8) luminosity; (9)x 2 s ■ 
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Fig. 6. Isodensity contours (logarithmically spaced) of the 55 clusters with significant structures. The title lists the coordinates of 
the center. The orientation is East to the left, North to the top. Galaxies belonging to the systems detected by DEDICA are shown as 
dots of different colors. Black, light green, blue, red, magenta, dark green are for the main system and the subsequent substructures 
ordered as in Table I A. ll Large symbols are for galaxies with My < -17.0 that lie where local densities are higher than the median 
local density of the structure the galaxy belongs to. Open symbols mark the positions of the first- and second-ranked cluster galaxies, 
BCG1 and BCG2 respectively. 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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Fig. 6. (continued) 
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